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Evolution of neutron star + He star binaries: an alternative 
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ABSTRACT 

It is difficult for intermediate-mass X-ray binaries to form compact intermediate-mass 
binary pulsars (IMBPs) with a short orbital-period (< 3 d), which have a heavy (> 
0.4 Mq) CO or ONeMg white dwarf companions. Since neutron star -I- He star binaries 
may experience common-envelope evolution, they have some advantage to account for 
the formation of short orbital-period IMBPs. In this work, we explore the probability 
of IMBPs formed by this evolutionary channel. Using Eggleton's stellar evolution code, 
considering that the dead pulsars were spun up by the accreting material and angular 
momentum from the He star companions, we have calculated the evolution of a large 
number of neutron star -|- He star binaries. Our simulated results indicate that, the 
NS -I- He star evolutionary channel can produce IMBPs with a WD of ^ 0.5 — 1.1 Mq 
and an orbital period of 0.03 — 20 d, in which pulsars have a spin-period of 1.4 — 200 
ms. Comparing the calculated results with the observational parameters (spin period 
and orbital period) of 9 compact IMBPs, the NS + He star evolutionary channel can 
account for the formation of 4 sources. Therefore, NS + He star binaries offer an 
alternative evolutionary channel to compact IMBPs. 
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1 INTRODUCTION 

Millisecond pulsars (MSPs) are characterized by weak sur- 
face mag netic fields {B ~ 10* G), short spin periods (1 < 
Pa < 30 ms), and nearly circular orbits with low-mass 
(0.15 Mq < Mwd < 0.4 Mq) He white dwarf (WD) compan- 
ions. These systems are generally called low-mass b inary pul- 
sars (LMBPs; e.g. iTauris fc van den Heuvelll200d ). In con- 
trast to LMBPs, observations fou nd an intermediate-mass 
binary pulsar (IMBP) po pulation l|Camilo et al.lll996l . [200ll : 
lEdwards fc Baileg|2001al ). which consists of a pulsar with a 
spin period of ~ 10 — 200 ms and a heavy (Mwd ^ 0.4 Mq) 
CO or ONeMg WD. 

At present, the formation mechanism of LMBPs is 
well understood. In the standard recycling scenario, LMBPs 
are the evolutionary products of low-mass X-ray binaries 
(LMXBs, with hydrogen-rich donor stars of < 1.5 A/0) via 
Cas e A or B Roche-lobe overflow (RLOF) (for a review , 
see iBhattacharva fc van den Heuvej Il99ll : IVerbuntI Il993h . 
A pulsar that passed through the so-called death line is 
spun up to miUisecond period by the accretion of material 
and angular momentum when the donor star fills its Roche - 
lobe l|Alpar et al.lll983 : lRadhakrishnan fc Srinivasan I1982D . 
Meanwhile, mass accretion onto a neutron star (NS) leads 
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to m agnetic field decay (|Taam fc van den Heuvell Il986l : 
iRoman i 1990), and the tidal interaction between the NS 
and the con vective envelope of the donor star induces a cir- 
cular orbit (|Phinnev 111993 : iPhinnev fc Kulkarni|[l99i) . The 
orbital-periods of LMBPs usually range from 0.2 days to 
several hundreds d ays, and ar e strongly correlated with the 
WD masses (.Tauris fc Savoniie 1999). 

Comparing with the LMBPs, IMBPs have relatively 
heavy WDs. In addition, the spin periods, the inferred sur- 
face magnetic fields, and the orbital eccentricities of the 
IMBPs are considerably higher than those of the LMBPs 
(L1.2003 ). These properties imply that the progenitors and 
the evolutionary history of IMBPs are distinct from those 
of LMBPs. It is generally thought that a large fraction of 
IMBPs are descendants of intermediate-mass X-ray bina- 
ries (IMXBs), in which the hydroge n-rich donor stars hav e 
a mass in the interval ~ 1.5 — 10 Mq (|van den Heuve]||l975l '). 
For example, IMBPs with an orbital-period of 3 — 50 d 
evolved from IMXBs with a donor star of ~ 2.5 — 5.0 Mq and 
an orbital-period of 3 — 10 d via Case B RLOF (Tauris et al.l 
I2OOOI : iPodsiadlowski et"ai]|2002l : pTauris et al.|[2012l ). In the 
IMXB evolutionary stage, the mass transfer rate is much 
larger than the Eddington accretion rate of the NS. After 
the rapid mass transfer, the donor star evolves into a low- 
mass star, and the binary begins LMXB evolutionary phase. 
The subsequent mass transfer rates are very low, and trigger 
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accretion disk instability. The accretion process shows short- 
lived outbursts separated by long-term quiescen ce, and the 
accretion efficiency markedly declines (|Lil I2OO2I ) . Although 
the NS in IMXBs only accretes a transferred amount of 
material of ~ 0.01 Mp) , it can still be sp un up to a short 
spm- period of ~ 10ms (iTauris et aklbOia ). 

The IMXBs evolutionary channel successfully explains 
the majority of IMBPs. However, once the IMXBs evolve 
into LMXBs phase, the binary orbits should widen be- 
cause the material is transferred from the less massive 
secondary to the more massive NS. Therefore, this evo- 
lutionary ch annel cannot prod uce short orbital-period (< 
3 d) IMBPs jTauris et al.ll2000l 'n. Nowadays, there exist 18 
known IMBPs (their orbital-periods range from 0.4 to 40 
d, see also Table 2), in which 9 compact IMBPs have an 
orbital-period of < 3 d. The origin of compact IMBPs still 
remains an unresolved puzzle. 

The CO WDs in IMBPs can be formed by tw o evo- 
lutionary channels as follows (|van den Heuve]|[T993 ): (i) A 
normal hydrogen-rich star burns helium and hydrogen shells 
around its degenerate CO core when it is on asymptotic gi- 
ant branch phase, (ii) A He star wi th a mass of < 2 Mq 
burns its helium shell l|Habets|[l983 ). The former route is 
an IMXB evolutionary channel, and the latter is a NS -I- He 
star evolutionary channel. Most of NS -I- He star binaries 
possess a compact orbit (their orbital-periods range from 
0.01 to 1.0 d, see also Figure 1 in Chen et al. 2011) because 
of the common-envelope evolution, and are potential pr ogen- 
itors of compact IMBPs. Recently, IChen et al.l (|201in sug- 
gested that the progenitor of the short orbital-period IMBP 
PSR J1802-2124 may be an NS -I- He star binary with a He 
star of 1.0 Ma and an orbital-period of 0.5 d. Meanwhile, 
iTauris et al.l (|2012l ) argued that, in post common envelope 
binaries case BB RLOF of He stars can form most binary 
millisecond pulsars with CO WD companions. In this work, 
we attempt to systemically investigate the initial parameter 
space of NS -I- He star binaries that could form IMBPs. 



2 DESCRIPTION OF BINARY EVOLUTION 
2.1 Stellar evolution code 

To exp lore the formation of dou ble NSs binaries, iDewi et al.l 
l|2002t ) and iDewi fc Poll (|2003l ) have calculated the evolu- 
tion of NS + intermediate-mass He star binaries. Similar to 
their work, we used an updated version of the stellar evolu- 
tion code developed by Eggleton (1971,1972,1973; see also 
Han et al. 1994; Pols et al. 1995) to calculate the evolution- 
ary sequences of NS (of mass Mns) + He star (of mass Mhc) 
binaries, until the systems evolve into the detached binaries. 
The initial chemical abundance of the He stars was taken to 
be Y = 0.98, Z — 0.02, and the ratio of the mixing length 
to the pressure scale height and the convective overshooting 
parameter were 2.0, and (Wang et al. 200^), respectively. 
In addition, the st ellar OPAL opacity table originated from 
a version given bv lChen fc 'fbiitl (|2007| ) (e. g. Rogers & Igle- 
sias 1992; Alexander & Ferguson 1994). 



^ As a possible path, an IMXB can evolve into a compact 
IMBP by anomalou s magnetic brak ing of Ap/Bp donor star 
djustham et al.ll2006l : IShao fc Li|[2012l) . 



Table 1. Three evolutionary phases of NSs during mass transfer. 



No. 


criteria 


evolutionary phase 


spin evolution 


1 


r-m < fc 


accretion phase 


spin-up 


2 


rm > rc 


propeller phase 


spin-down 


3 


rm > ric 


radio phase 


spin-down 



Note: rm,rc, and ric represent the magnetosphere radius, the 
co-rotation radius, and the light cylinder radius of the NS, 
respectively. 



2.2 Mass transfer and angular momentum loss 

During the evolution of binaries, orbital angular momentum 
loss play a vital role. In calculation, we consider three types 
of angular momentum loss from the binary system, which 
are described as follows: 

1. Gravitational wave radiation (GWR). For He star-NS 
binaries with a compact orbit, GWR can effectively carry 
away orbital angular momentum and leads to mass transfer. 
This angular momentum loss rate is 



^GR = 



32G^/2 M^sM^M 



1/2 



5c5 



aV2 



(1) 



where G and c are the gravitational constant and the speed 
of light, respectively; M — A/ns + Muc is the total mass of 
the binary, and a is the binary separation. 

2. Mass loss. During a rapid mass transfer, there may 
exist mass and angular momentum loss due to the limita- 
tion of Eddington accretion rate, in which MEdd — 3.0 x 
10~* M0yr~^ for He-rich accretion material. When the mass 
transfer rate | Mhb |> M^dd, we assumed that the mass 
loss rate of the system M = Mne + iVfcdd proceed s accord- 
ing to isotropic re-emission (^ Soberman et al.lll997l ). and car- 
ries away the specific orbital angular momentum of the NS. 
The orbital angular momentum loss rate by the isotropic 
re-emission is given by 



Jm 



M2 



a ft, 



(2) 



where Q is the orbital angular velocity of the binary. 

3. Magnetic braking. Even if there exist no mass ex- 
change, the coupling between the magnetic field and the 
stellar winds can also indir ectly carry away the orbit al angu- 
lar momentum of binaries l|Verbunt fc Zwaanlll98ll ). In this 
work, the induced magnetic braking model given by Sills et 
al. (2000) was adopted. 



2.3 Spin evolution of the NS 

In the stellar evolution code, the evolution of the spin and 
the surface magnetic field of NSs were also considered. In 
summary, we consider three spin evolution stages in Table 
1. The detail ed input physics f or this subsection see also 



i. ine detail ed input pnysics t c 
section 2.3 in lLiu fc Chgiil (120111 ) 



3 SIMULATED RESULTS 

To explore the initial parameter space of NS -I- He star bi- 
naries that can form IMBPs with various spin-periods, we 
have calculated the evolution of a large number of NS + 
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Age (Myr) 



Figure 1. Evolutionary tracks of an NS + He star binary with 



Mi 



I.OM0, and P^j. 



2.5 d. The solid and dashed curves 



represent the evolution of the NS mass and the He star mass in 
the top panel, the mass transfer rate and the orbital period in the 
middle panel, and the spin period and the surface magnetic field 
of the NS in the bottom panel, respectively. 
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Distribution of the initial He star mass and the 

of NS + He star binaries that can evolve 



Figure 2 

initial orbital period 
into IMBPs via the He stars evolutionary channel. The filled cir- 
cles, open circles, and filled triangles denote IMBPs with spin- 
period of 1 — 10 ms, 10 — 200 ms, and 200 — 1000 ms, respectively. 
The open squares correspond to NS + He star binaries that expe- 
rience dynamically unstable mass transfer, and the crosses repre- 
sent the binaries in which the He star does not fill its Roche lobe 
in the Hubble time. 



He star binaries. Employi ng the population synthesis ap- 
proach, IChen at al.l (2Q1^ found that the initial He star 
masses and the initial orbital-periods of NS -|- He star bina- 
ries are M|i<, ~ 0.5 - 2.0 Mq, and P^^b ^ 0.01 - 10 d, re- 
spectively. According to their results, we take the simulated 
grids to be M|i„ = 0.5 - 2.0 Mq, and P^^ = 0.06 - 10.0 d. 
For normal radio p ulsars, the maximum spin period is 11 s 
l|Manchesterl 1200^ 1 . and the final spin periods of NSs is in- 
sensiti ve to the initial sp in-periods and the initial magnetic 
fields (jWang et al.ll20l"ll l. Therefore, the initial spin-period 
and the initial magnetic field of the NSs are taken to be 10 s, 
and lO'^^ G, respectively. Furthermore, a canonical NS mass 
Mj^s = 1.4 Mq was adopted. 

In Figure 1, we show the detailed evolutionary se- 
quences of a NS -I- He star binary with an initial mass = 
1.0 M0, and an initial orbital period Pq^i, = 2.5 d. After the 
exhaustion of central He, the He star begins case BB mass 
transfer at the age of t ~ 18.0 Myr. The rapid transfer of 
He-rich material occurs at a high rate of ^ 10~® M0yr~^. 
Because of the Eddington accretion rate, a majority of the 
transferred material is ejected in the vicinity of the NS, and 
form the isotropic re-emission. When t — 18.16 Myr, the 
mass transfer ceases, and the binary becomes a detached 
system. Hereafter, the He star evolves into a CO WD and 
begins cooling. After mass transfer of 0.16 Myr, the accreted 
mass for the NS AMns = MEdd x 0.16 Myr = 0.0048 Mq. 
Though the accreted mass is ~ 3% of the transferred matter, 
it is sufficient to spin the NS up to Ps ~ 45 ms, and causes 
the magnetic field to decay to 2 x 10^" G. Owing to the less 
massive He star transfers the material to the more massive 
NS, the orbital period continuously increases to 3.53 d. 

Figure 2 compiles the final evolutionary fate of NS + 
He star binaries in the initial He star mass - orbital period 
(M^g — Pojb) plane. It is seen that NS -I- He star binaries with 
He stars of 0.8 - 1.4 Mq can evolve into IMBPs. When the 
initial orbital period Pl-^^, ^ 0.1 d, most evolutionary prod- 



ucts are IMBPs with a millisecond period. If Pl^-^a > 0.1 d, 
nearly all NSs are only mildly recycled, and their spin-period 
are in the range of 10 -200 ms. In particular, in our calcu- 
lated grids, there exist 4 NS + He star binaries that pro- 
duce slow spin pulsars with a spin period of 200 — 1000 ms. 
When the mass ratio of binaries q = Mhc/Mns > 1, the 
envelopes of the He stars is dominated by convection at the 
onset of RLOF (Ivan den HeuvelllioogI) . Therefore, a run- 
away mass transfer is triggered, and nearly all NS -I- He star 
binaries with q > \ would experience CE evolution. How- 
ever, for H main sequence companions in LMXBs or IMXBs 
with radiative (or only slightly convective) envelopes, and 
the mass transfer is stable, and can avoid CE evolution- 
ary stage even if </ > 1 (see also Tauris fc Sayonii j 19991: 
Tauris et all l200d: IPodsiadlowski et al.l I2OO2I : iLiu fc ChenI 



201ll : IShao fc Lill2012l ). 



In order to study how the accreted mass influences the 
spin evolution of NSs, in Figure 3 we plot the evolution- 
ary results in the accreted mass vs. the final spin period of 
NSs (AMns — Ps) diagram. It is clear that, if NSs accrete 
He-rich material of > 0.02 Mq , they will be spun up to mil- 
lisecond period, and ev olve into IMBPs like PSR J1614-2230 
dPemorest et al.ll2010l l. As shown in this figure, there exist 
a practically linear logarithmic inverse correlation between 
the final spin-period and the accreted mass of NSs, which 
originates from our model assumption for the spin-up of NSs 
(see equation 9 in Liu & Chen 2011). 

In Figure 4,we present the distribution of our simulated 
results in the Porb,f — Ps diagram. One can see that our NS 
-I- He star evolutionary channel can form IMBPs with an 
orbital period of 0.03 - 20 d and a spin-period of 1.4 — 1000 
ms. It is difficult for our evolutionary scenario to produce 
IMBPs with a relatively long orbital-period (> 20 d) , which 
should be derived from the selection of the distribution of 
the initial orbital-periods. In Table 2, we list the observed 
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Table 2. Observed parameters for 18 IMBPs. 



Pulsars 


P^(ms) 

J. g yiiioy 


-'orb y^'^j^J 


References 


J1420-5625 


40.3 


34.1 


1 


J1810-2005 


32.8 


15.01 


2 


J1904+0412 


71.1 


14.93 


2 


J1454-5846 


45.2 


12.42 


2 


J1614-2230 


3.15 


8.69 


3,4 


J0621+1002 


28.9 


8.319 


5,6 


J1022+1001 


16.5 


7.805 


6 


J2145— U75U 


16.1 


6.839 


7 


J1603-7202 


14.8 


6.309 


8 


J1157-5112 


43.6 


3.507 


9 


J1528-3146 


60.8 


3.18 


10 


J1439-5501 


28.6 


2.12 


11,12 


J1232-6501 


88.3 


1.86 


2,13 


J1435-6100 


9.35 


1.35 


2,13 


B0655+64 


195.7 


1.03 


14,15 


J1802-2124 


12.6 


0.699 


11,16 


J1757-5322 


8.87 


0.453 


17 


J1952+2630'= 


20.7 


0.392 


18 



Refere nces. QlHo bbs et aTl ll2004h : f2lCamilo et al 
(2001); ('3 lCrawford et al l bOOdt) -. f4 lDemorest et al, 
f5lSplaver et al.l tOO^: f6lCarnilo et al, 
(7lBailes et al.l lll994h : C8 lLorimer et al.l lll99l " 
esl tOOlJ ): d Oljacobv et al.l 



("9' Edward s &: Baile: 
I'll' Faulkner et al.l 
(13' 

(is: 
(17: 



Manchester et a' 




Edwards fc Baile j l|2001bh : (18 iKnisDel et al.l pOllF 



(12 lLorimer et al-. 
(1 4ljones fc Lvn3 
(161Ferdmar^_et_^ 
1 



parameters for 18 IMBPs. To compare with observations, 
we also show the locations of 18 IMBPs by the filled stars in 
Figure 4. As shown in this figure, our calculated results can 
roughly fit 5 - 6 observed sources. Other sources have much 
shorter spin period than that predicted by the model. How- 
ever, for 9 short orbital-period IMBPs, our evolution results 
can fit 4 sources, suggesting that the NS + He star evolu- 
tionary channel may be an interesting evolutionary path to 
short orbital-period IMBPs. 



0.1 - 



0.01 - 



lE-3 - 



lE-4 




Figure 3. Distribution of our simulated results in the accreted 
mass of NSs vs. the final spin-period (AM^s — Ps) diagram. 



100 




0.1 r 



Figure 4. Distribution of our simulated IMBPs in the final or- 
bital periods vs. the spin-period of NSs (Porb.f ~ Ps) diagram. 
The solid stars represent the observed sources. 



4 DISCUSSION AND SUMMARY 

In this paper, we attempt to investigate if IMBPs can have 
evolved from NS -f He star evolutionary channel. We have 
performed numerical calculations for the evolution of NS -I- 
He star binaries consisting of a 1.4 Mq NS and a 0.5 - 2.0 Mq 
He star companion. Our main results can be summarized as 
follows. 

1. When the He star mass is in the range of 0.8 — 1.4 Mq, 
NS -I- He star binaries can evolve into IMBPs. If the initial 
orbital period P^^t, is too short (< 0.1 d), the NS will be fully 
recycled, and can be spun up to < 10 ms. In other cases, 
a mildly recycled pulsar with a spin-period of 10 -1000 ms 
would be produced. 

2. During the mass exchange, the accreted mass of most 
NSs is in the range of 0.001 — 0.18 Mq, which is sufficient to 
spin up the NS to a period of 1 — 200 ms. An accreted mass 
of 0.02 Mq is the threshold that the NSs are recycled to be 
millisecond period. 

3. The NS + He star evolutionary channel can pro- 



duce IMBPs with an orbital-period of 0.03 - 20 d. Out of 9 
known short orbital-period IMBPs, our model can address 
the formation of 4 compact sources. Comparing with the 
IMXBs evolutionary channel, the NS -I- He star binaries can 
form short orbital-period (< 3 d) IMBPs wi thout invoking 
anomalous magnetic b raking of Ap/Bp stars jjustham et al.l 
l2006l : IShao fc Lill2012[ l 

4. Comparing with the observations, our simulated 
IMBPs tend to have a long spin-period, or a short orbital- 
period. There may be three possibilities to explain the dif- 
ference between observations and our model predictions. 
Firstly, the NS -I- He star evolutionary channel only pro- 
duces a part of obser ved IMBPs, and t he remainder de - 
scended from IMXBs (|Tauris et al.ll2000l : IShao fc Lill2012l '). 
Secondly, the theoretical model of magnetic braking needs 
to b e revised, or sup er-Eddington accretion of NSs may oc- 
cur l|Begelmanl 120021 '). Thirdly, in our calculation we take a 
canonical NS mass (of mass 1.4 Mq), while the NS masses 
may have a large range of values. Actually, the initial pa- 
rameter space (especially initial orbital-period) for form- 
ing IMBPs in Porb — diagram increases with the in- 
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crease of the NS mass (ishao fc Lil ioii). Observational 

and theoretical studies indicate that NSs formed by iron- 
core collapse supernovae and electron-capt ure supernovae 



Timmes et alj 199e 


2005 


: 1 Schwab et al 


2012 


). For example 



Heger et all |2003|: 



I2OIOI : IZhang et al.l I2OIII : lOzel et al 



Wooslev fc Janka 



Isar PSR J1614- 2230 
with a mass of 2.0 Mq was proposed to be originated from 
IMXB with a massive NS of 1.6 - 1.7 Mq l|Lin et al.ll2"oill : 



iTauris, Langer fc Krameill201ll '). 
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